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Measurements of resistivity and impurity concentration in heavily doped 
silicon are reported. These and previously published data are incorporated 
in a graph showing the resistivity (at T = 300°K) of n- and p-type silicon 
as a function of donor or acceptor concentration. 

The relationship between surface concentration and average conductivity 
of diffused layers in silicon has been calculated for Gaussian and comple- 
mentary error function distributions. The residts are shown graphically. 
Similar calculations for subsurface layers, such as a transistor base region, 
are also given. 

I. INTRODUCTION 

A diffused layer in silicon is generally characterized by four parame- 
ters: the concentration, C, , of diffused donors or acceptors at the surface, 
the concentration, C a , of acceptors or donors originally in the material 
(background concentration), the depth, xj , of the resultant junction, 
and the sheet resistivity, p, , of the layer. A knowledge of the relationship 
between these parameters is essential to the establishment of device 
processing recipes, the evaluation of diffusion techniques, and investiga- 
tions of the thermodynamic properties of silicon. 

The desired relationship may be readily calculated, given a knowledge 
of the distribution of the diffused impurities, the variation of the re- 
sistivity of n- and p-type silicon with donor or acceptor density, and a 
fast electronic computer. The results of such a computation were first 
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made generally available three years ago, in the form of curves relating 
C, to 1/psXj for a given C B , for n- and for p-type layers in silicon, and 
for several common distributions. 1 Recent calculations, however, based 
on new and more extensive silicon resistivity data, have indicated con- 
siderable error in the earlier results. Thus a comprehensive recomputa- 
tion has been undertaken, the outcome of which is presented herewith. 
A necessary adjunct to the calculation is an accurate knowledge of the 
resistivity of n- and p-type silicon with varying dopant concentration. 
To this end, most of the extant data have been reviewed and supplemented 
here and there with some new determinations. The results of this search 
are also presented here. 

II. THE RESISTIVITY OF SILICON AS A FUNCTION OF IMPURITY CONCEN- 
TRATION 

The variation of the resistivity of silicon at 300°K as a function of the 
concentration of acceptors or donors is shown in Fig. 1. This graph 
represents the author's judgment of a most reasonable compromise to 




IMPURITY CONCENTRATION (CM -1 ) 

Fig. 1 — Resistivity of silicon at 300°K as a function of acceptor or donor con- 
centration. 
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Table I — Resistivities and Impurity Concentrations 
in Silicon (T = 300°K) 



Resistivity 


Impurity 


Impurity 


Carrier 


(ohm-cm) 


Concentration (cm -3 ) 


Concentration (cm -3 ) 


0.00076 


B 


1.66 X 10 20 




0.00089 


B 


1.41 X 10 20 




0.0010 


B 




1.49 X 10 20 


0.0010 


B 


1.12 X 10 20 




0.0012 


B 


1.04 X 10 20 




0.0011 


B 


1.12 X 10 20 




0.0014 


B 


9.23 X 10 19 




0.0013 


B 


8.84 X 10 19 




0.0007 


B 


1.43 X 10 19 




0.0073 


B 


1.43 X 10 19 




0.013 


B 


7.41 X 10 18 




0.014 


B 


7.03 X 10 18 


■ 


0.00095 


As 


1.80 X 10 20 




0.00094 


As 


1.86 X 10 20 




0.00094 


As 




1.1 X 10 20 


0.00093 


As 


1.87 X 10 20 




0.00094 


As 


1.97 X 10 20 




0.00088 


As 


2.10 X 10 20 




0.00088 


As 


2.19 X 10 20 




0.00089 


As 




1.1 X 10 20 


0.00083 


As 


2.30 X 10 20 




0.00083 


As 


2.20 X 10 20 




0.00080 


As 


2.46 X 10 20 




0.00082 


As 


2.44 X 10 20 





the mass of available and not altogether compatible data on the subject. 
These data include most of the previously published work (Refs. 3-12), 
recent, unpublished results kindly provided by other investigators, 2 ' 13 
as well as some measurements obtained expressly for the present study. 

The last data are shown in Table I. The crystals involved were pulled 
from quartz crucibles, and hence can not be expected to be particularly 
low in oxygen content. After dissolution of the boron-doped crystals 
and separation of the dopant, 1 boron concentrations were determined by 
a photometric carmine technique essentially similar to published meth- 
ods. 1 Arsenic concentrations were measured by gamma-ray spectrometry 
after pile neutron activation. Resistivity measurements were done with 
a four-point probe. In the case of a few samples, resistivity and carrier 
concentration were measured in Hall-effect apparatus (where it was 
assumed mw/m = 1 ) • 

Drawing curves through these many points was accomplished by a 
succession of smoothing procedures, which were primarily visual. 75 per 
cent of the data points deviate less than 10 per cent from the curves thus 
obtained, both for the p-type and the n-type cases. The uncertainty is 
greatest in the degenerate region. For p-type silicon, suitable data be- 
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come scarce at dopings greater than 10 19 cm -3 , and none are available 
beyond 3 X 10 20 cm -3 . For n-type material, there is an abundance of 
rather conflicting data representing donor concentrations between 10 Ifl 
cm" 3 and 6 X 10 20 cm -3 . In this region a 10 per cent variation in the 
chosen line still includes 67 per cent of the data, however. 

A single pair of curves obviously can not characterize with the same 
degree of accuracy all silicon material, regardless of dopant employed 
or degree of compensation. However, over the range 10 cm - ^ Ni ^ 
10 20 cm -3 , and subject to the limitations discussed below, Fig. 1 is con- 
sidered to be within 10 per cent of reality. This graph refers specifically 
to uncompensated silicon containing a donor or acceptor impurity con- 
centration, Ni , consisting of arsenic, phosphorus, or antimony for 
n-type, and aluminum, boron, or gallium for p-type material. (Actually, 
even among samples doped with the aforementioned impurities, small 
but consistent differences in carrier concentration and mobility, depend- 
ing on the specific choice of donor or of acceptor, have been reported 
recently for silicon in the 0.001 ohm-cm region. 10,12 ) In case of moderate 
compensation, the net impurity density, \N A — N D \, should be used 
for Ni . However, heavy compensation requires allowance for the added 
impurity scattering. 

For impurity densities near or greater than 10 2 cm - ', Fig. 1 can not 
be considered very reliable. At such concentrations, impurity band 
conduction is prominent and its effects are apt to differ appreciably 
depending on choice of impurity. Even more serious are the degrees of 
impurity precipitation and lattice imperfection which occur in highly 
doped material and which furthermore vary with growth conditions 
and history of the crystal. It will be noted with some consternation that 
the p-type and n-type curves are shown to cross near Ni = 3 X 10 cm . 
The paucity of data, of course, casts considerable doubt on this result. 
However, for what they are worth, such are the indications. Perhaps this 
can be understood in light of the acceptor action of imperfections, 
especially vacancies, which are abundant in very highly doped material. 

The calculations discussed in the remainder of this paper require a 
mathematical representation of Fig. 1. Straight-line approximations of 
the form (1/p) = BNi" have been obtained, which depart 10 per cent 
from the desired curve at the turning points and rapidly approach 
coincidence elsewhere. The parameters B and a are listed in Table II 
for the respective straight-line regions. 

III. DIFFUSION PROFILES AND CALCULATIONS 

The diffusion profiles of current practical interest are the comple- 
mentary error function, C x = C„ erfc (x/2\/Dt), and the Gaussian, 
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Table II — Values of B and a in the Equation (1/p) = BNi a , 

Representing Straight- Line Approximations to the p vs 

Ni Curves of ii-type and p-type Silicon (T = 300°K) 



Region (cm 3 ) 


B 


" 


n-lype 








2.35 X 10 20 


^ N D 


1.04 X 10-° 


0.456 


6.00 X 10 19 


^ N D ^ 2.35 X 10 2 ° 


1.43 X 10~ 12 


0.744 


9.50 X 10 18 


< N D < 6.00 X 10 19 


2.00 X 10- 18 


0.940 


1.00 X 10 17 


<, N D ^ 9.50 X 10 18 


6.93 X 10" 9 


0.543 


3.50 X 10 15 


<, N D < 1.00 X 10 17 


6.97 X 10~ 14 


0.837 




N D g 3.50 X 10 15 


2.00 X 10- 16 


1.000 


p-type 








1.50 X 10 1D 


^ N A 


4.00 X 10- 17 


0.966 


2.40 X 10 18 


< N A < 1.50 X 10 la 


1.47 X 10 -14 


0.832 


1.50 X 10 16 


< N A < 2.40 X 10 18 


3.30 X 10" 11 


0.650 




N A ^ 1.50 X 10 1 " 


7.20 X 10" 17 


1.000 



C x = C s exp ( —x 2 /4Dt). In these expressions, x, D, and t are the depth, 
diffusion coefficient (assumed independent of impurity density), and 
time, respectively. C x is the concentration of the diffused impurity at 
depth x and C 8 , that at the surface. The former distribution is expected 
when diffusion takes place with the surface concentration C 8 held con- 
stant; the latter when the total impurity diffusing is constant. Unfor-. 
tunately it must be admitted that the accuracy of these expectations is 
open to question in some situations."' 1 Also, precipitation and compen- 
sation of impurities near the surface may further distort the distribution. 
However, it is still useful to solve the problem under these assumptions, 
leaving corrections for later determination. 

The "average conductivity" of a diffused layer (which throughout 
this paper is assumed to be diffused into a silicon slice of opposite con- 
ductivity type and uniform doping C B ) is given by the expression 



a = 1/psXj = (l/xj) f ' qnC 
Jn 



dx 



where q is electronic charge, ju the carrier mobility typical of a total 
ionized impurity density of C x + C B , C = r(C x — C B ) is the density of 
carriers, r being the fraction of uncompensated diffused impurity atoms 
which are ionized, and C x the total density of diffused impurity atoms at 
depth x. (Possible variation of the mobility as aJ:uiic±ion of the proxim- 
ity of the surface is a hazard which should be recognized in passing but 
is otherwise ignored in the present calculation. ) Multiplying and dividing 
within the integrand by r' (C x + C B ), where r' is the ionized fraction 
associated with an uncompensated dopant density of (C x -+- C B ), and 
writing 
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qtir'(C x + C B ) = a(cx+cB) = B(C X + C B ) a 
the average conductivity becomes 

a = Ufa) f' (r/r')(C x - C B )B(C X + C^dx. 

Now (r/r') represents the ratio of degrees of ionization corresponding 
to C x — C B and C x + C B respectively. This ratio is very nearly unity 
unless C x and C B are comparable in magnitude. Such is the case only for 
the lamina nearest the junction, which contributes negligibly to the 
conductance of the whole layer. Hence, {r/r') may be justifiably taken 
as equal to unity, and writing C x = C s f(x), where /(.t) depends on the 
profile of interest, 

a = (l/xj) [*' [C s f(x) - C B ]B[C s f(x) + C B r l dx. 
Jo 

A program for the evaluation of this expression has been devised 
previously by others and employed in the analysis of diffused layers in 
germanium. 17 With slight additions to facilitate automatic plotting, the 
same program has been used in the present work. Computations were 
performed on an IBM 704, and plotting of points was carried out with 
an Electronic Associates Variplotter. 

IV. PRESENTATION OF RESULTS 

Of frequent interest in transistor design and in the analysis of diffused 

layers, are the characteristics of a "subsurface" layer such as illustrated 

in Fig. 2. This layer, bounded on one side by the junction and on the 

other by a plane paralleling the junction at depth x, may be characterized 

by an average conductivity 

1 f x ' 

5 = l/W(xj - x)] = / qtiC dx 

\Xj — X) J x 

where p/ is the sheet resistance of the subsurface layer. It will be recog- 
nized that the base region of a diffused-base, alloyed-cmitter transistor 
is an example of a subsurface layer. Another example is that portion of 
a diffused layer remaining after removing the top strata of depth x. 
Here, however, it must be remembered that the value of C a specifying 
this layer pertains to the original surface at x = 0. 

Since a subsurface layer becomes the entire diffused layer when x = 0, 
it is convenient to display the properties of both in the same plot by 
introducing the parameter (x/xj). On pages 394 to 410 such graphs are 
presented for n- and p-type diffused layers of Gaussian and comple- 
mentary error function profile. Each graph contains the family of ten 
curves (x/xj) = 0, 0.1,- • -,0.9, and relates the average conductivity of 
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DEPTH — > 
Fig. 2 — Profile of a diffused layer with subsurface layer shaded. 

each layer to the surface concentration (at the original surface) for a 
given value of C B . A separate graph is required for each value of C B , 
which in the present work ranges from 10' cm - ' to 10" cm - ' at one- 



decade intervals. In each 

■»21 



plot the range of surface concentrations 

-3 



spanned is from C B to 10" cm - '. The so-called "Backenstoss" curve for 
a particular C B is simply the right-most line (x/xj = 0) in each graph. 
The wiggle in the n-type average conductivity for diffusant concentra- 
tions near 10 19 cm" is ascribable to the rather large change in slope oc- 
curring in the n-type resistivity plot at Ni = 10 1 cm - . 
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Fig. 3 — Average conductivity of n-type complementary error function layers 
silicon. 
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Fig. 3 (cont.) — Average conductivity of n-type complementary error function 
layers in silicon. 
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Fig. 3 (cont.) — Average conductivity of n-type complementary error function 
layers in silicon. 
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Fig. 3 (cont.) — Average conductivity of n-type complementary error function 
layers in silicon. 
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Fig. 3 (cont.) — Average conductivity of n-type complementary error function 
layers in silicon. 
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Fig. 3 (cont.) — Average conductivity of n-type complementary error function 
layers in silicon. 

397 





[X/Xj)n 


0.9 






0.8 




0.7 


0.6 


o.! 




0.4 0.3 


0.1 
























/ 






/ 






A ' 


f //// 






























/ 








— 7 




/ //// 






















/ 






/ 






/ 




























/ 






/ 




































/ 












































































































C B = I0 20 CM- 3 
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AVERAGE CONDUCTIVITY, ff = l/[/5' s (Xj-X)] (OHM-CM)" 1 

Fig. 3 (cont.) — Average conductivity of n-type complementary error function 
layers in silicon. 
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Fig. 4 — Average conductivity of n-type Gaussian layers in silicon. 
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Fig. 4 (cont.) — Average conductivity of n-type Gaussian layers in silicon. 
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Fig. 4 (cont.) — Average conductivity of n-type Gaussian layers in silicon. 
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Fig. 4 (cont.) — Average conductivity of n-type Gaussian layers in silicon. 
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Fig. 4 (cont.) — Average conductivity of n-type Gaussian layers in silicon. 
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Fig. 4 (cont.) — Average conductivity of n-type Gaussian layers in silicon. 

401 



I0 2 ' 


(X/Xj) = 


0.9 








0.8 







7 


0.6 0.5 0.4 


0.2 


























/ 


' 




f- 




/ 




t 


-rrt 


rtr 
























/ 








/ 




/ 




/ 


U- / / 


// 
























/ 








P 




/ 




' 






























/ 
















/ 


l ///// 








































/ 


' 


'///y 




















































C B = l02"CM-3 




















































/ 


' /// 






































/ 


/ 


^ 


? 


V 
































^ 


£ 


^ 


^ 


^ 


V 
































\0'° 















































* i I 



4 s e 



AVERAGE CONDUCTIVITY, 0° = l/[/>' s (Xj-X)] (OHM-CM)"' 

Fig. 4 (cont.) — Average conductivity of n-type Gaussian layers in silicon. 

(X/Xj) = 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 
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Fig. 5 — Average conductivity of p-typc complementary error function layers 
in silicon. 
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Fig. 5 (eont.) — Average oonduetivity of p-type complementary error function 
layers in silicon. 
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Fig. 5 (cont.) — Average conductivity of p-type complementary error function 
layers in silicon. 
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Fig. 5 (cont.) — Average conductivity of p-type complementary error function 
layers in silicon. 
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Fig. 5 (cont.) — Average conductivity of p-type complementary error function 
layers in silicon. 
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Fig. 5 (cont.) — Average conductivity of p-type complementary error function 
layers in silicon. 
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AVERAGE CONDUCTIVITY, <T = l/[/J^ (Xj-X)J (OHM-CM)" 1 

Fig. 5 (coot.) — Average conductivity of p-t-ype complementary error function 
layers in silicon. 
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Fig. 6 — Average conductivity of p-type Gaussian layers in silicon. 

406 



RESISTIVITY IX SILICON 



407 



(X/rj) = 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 




2 3 4 1 ea |0 -, J 3 456. , 2 3 4 5 6 e |£) , 2 3 4 5 «. |q2 2 3 4 56.^ 2 J 4 S I « |tf( 

AVERAGE CONDUCTIVITY, ff=\/[p' s (Xj-X)] (OHM-CM)"' 



Fig. 6 (cont.) — Average conductivity of p-type Gaussian layers in silicon. 
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Fig. G (cont.) — Average conductivity of p-type Gaussian layers in silicon. 
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Fig. G (cont.) — Average conductivity of p-type Gaussian layers in silicon. 
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Fig. 6 (cont.) — Average conductivity of p-type Gaussian layers in silicon. 
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Fig. G (cont.) — Average conductivity of p-type Gaussian layers in silicon. 
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Fig. 6 (cont.) — Average conductivity of p-type Gaussian layers in silicon. 



